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Random  Covering  and  Packing  on  the  Line 

by  Howard  J.  Veiner 

University  of  California,  Davis  and  Stanford  University 


I.  Introduction:  Covering  Model.  Unit  length  segments  are  sequentially 
placed  independently  at  random  uniformly  on  a  line  segment  [0,x)  until 
that  segment  is  completely  covered.  The  covering  models  of  Renyi  and 
Solomon  will  be  considered  in  section  2,  and  asymptotic  results  for  the 
two  dimensional  cases  for  both  models  will  be  given  in  section  3.  The 
variance  and  other  extensions  will  be  indicated. 

A  packing  and  covering  model  for  the  one  dimensional  abacus  grid 
will  be  indicated.  The  one -dimensional  Renyi  model  is  as  follows. 

The  center  of  each  unit  segment  is  chosen  uniformly  on  [0,x).  The 
unit  segment  so  chosen  is  placed  on  [0,x),  x  >  1.  The  process  is  re¬ 
peated,  and  a  new  unit  segment  is  kept  (parked)  if  and  only  if  its 
center  is  on  an  uncovered  portion  of  [0,x),  otherwise  it  is  discarded. 

The  process  stops  as  soon  as  the  [0,x)  is  completely  covered.  All 
unit  segments  except  those  covering  the  ends  must  lie  within  [0,x>. 
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Let  M(x)  “  mean  o£  the  minimum  number  of  unit 
segments  covering  [0,x)  in  the  Renyi  model. 


(1.1) 


The  Solomon  model  allows  the  center  of  each  unit  segment  to  be  chosen 
uniformly  on  [0,x).  Unit  segments  are  placed  without  exclusion  or  further 
constraint  until  the  [0,x)  is  covered.  Let 


R(x)  *■  mean  of  the  minimum  number  of  unit  segments 
covering  (0,x)  in  the  Solomon  model. 


(1.2) 


Results  for  a  model  related  to  the  Solomon  model  (1.2)  are  obtained 
in  Solomon  (1966).  The  results  here  are  by  different  considerations. 

The  parking  models  (Blaisdell  et  al  (1970),  Solomon  (1966)  are 
related  to  these  problems  and  references  relevant  to  sections  I,  II, 

VII  are  contained  there. 

let,  in  either  model. 


X^(t),  ^(t)  be  I.I.D.,  each  distributed  as  X(t), 

which  is  the  minimum  total  number  of  unit  length  (1.3) 

cars  needed  to  cover  [0,t). 


} 

Then  conditional  on  the  first  placement  with  center  at  t  +  — ,  one 
may  represent  X(x)  as 


X(x)  -  Xx(t)  +  X^x-t-l)  +  1  (1.4) 

Taking  expectations  of  (1.4)  for  the  Renyi  model  and  combining 


two  Integrals  on  the  right,  one  obtains 
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X-i 

M(x)  -  1  +  §  f  2M(u)du,  x  >  j 
X  v0  ,  2 
and  M(x)  ■  1,  0  <  x  <  j • 


(1.5) 


Hence  for  x  >  1, 


X_1 

M(x)  -  1  +  J  +  |  i  2  M(u)du. 


(1.6) 


Lemma  1.  If  A(x),  x  >  1,  satisfies 

1 

4<x)(<)  I  (  2l(u)du 


(1.7) 


and 


then 


A(x)(|)  0,  for 
1  <  xQ  <  x  <  xQ  + 


M*)  (~)0  for  all  x  >  xQ. 


(1.8) 


(1.9) 


Proof.  Assume  the  upper  inequality  in  (1.7)  and  (1.8).  Then 
(1.8)  into  (1.7)  for  Xq  +  j<x<Xq  +  1  yiel^s 

A(x)  -  \  [J*°  A(u)du  +  jX  2  A(u)du] 
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(1.10) 


The  first  term  on  the  right  of  (1.10)  is  non-negative  since 
A(x0  +  2 )  >  0  and  the  integrand  of  the  second  term  is  non-negative  by 
(1.8).  The  result  now  holds  for  x^  <  x  <  Xq  +  1,  so  that  the  process  is 
continued  as  before,  and  an  induction  completes  the  proof.  The  lower 
inequality  is  similar. 


Theorem  1.  For  x  >  1, 


1.485  x  +  .485  <  M(x)  <  1.5x  +  .5 


Proof.  An  explicit  computation  verifies  that  if 


Lt(x)  b  1.485x  +  .485 


L j (x)  a  1 , 5x  +  .5 


then  for  x  >  1, 


i.2  rx“2 


L.  (x)  <  1  +  -  +  -  j  L.  (u)du 


x  xvl  “l' 
2 

xJ. 

L2(x)"1  +  x  +  x^ 

2 


and  for  1  <  x  <  j, 

Lj  (x)  <  M(x)  <  L2  (x)  . 
Hence,  defining,  for  i«l,2, 

\(.x)  -  L^x)  -  M(x), 

it  follows  that 


K-I 


*!<*>  sLr  '\M4» 


1 

2 


V«)  -  f  j*  \(u)du 


(1.11) 

(1.12) 

(1.13) (i) 

(1.13)  (li) 

(1.14) 

(1.15) 

(1.16) (i) 

(1.16)  (ii) 
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and  for  1  <  x  <  j,  it  is  easily  computed  that 


(1.17) 

^(x)  >  0. 

The  theorem  now  follows  from  (1.14)  -  (1.17)  and  the  lemma. 

Taking  expectations  of  (1.4)  for  the  Solomon  model,  one  obtains 
for  x  >  0, 

R(X>  "  1  +  jin  j^R<u)du  (1.18) 

with  R(0)  -  1. 

This  Volterra-type  equation  has  a  unique  solution,  which  is  easily 
seen  to  be,  for  x  >  0 

R(x)  ■  2x  +  1.  (1.19) 


To  obtain  precise  asymptotic  behavior  of  M(x) ,  define  the  Laplace 
transform  for  0  <s  <  1 

L(s)  s  J  e”sxM(x)dx.  (1.20) 

2 

Since  M(x)  <  2x  +  1,  L(s)  is  well  defined.  From  (1.5)  it  immediately 
follows  that 

/  .  />i  i  o'. 

(1.21) 


%  .  2e  %  2e  -s/2  /'l  .  3 

L  <8>  + 1 —  L(8)  -  —  - e  j;  + 1 ) 

s  s 


with  solution,  using  the  fact  that  M(x)  <  cx. 


L(s) 


f*e"u/2  f 

'  *  “  J- 


c"  -u/2 

2J 

v 


r  -v/2  '  1  .  3  N  2e“V  |  .  n 

L  e  \2t  +  -2  >T  J  (l 


.22) 
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Since  for  0  <  s  <  1, 

-u/2 


J 


n  8/2  -U  - 

du  -  -  Y  -  J  ~~  du  -Wb, 

o 


(1.23) 


for  0  <  s  «  1,  it  is  seen  by  a  Taylor  expansion  to  one  term. 


T  /  \  ~  c  ,  c  -  1 
L(s)  "  “J  +  — 

s 


(1.24) 


where 


(1.25) 


lim  x-1M(x)  »  c,  (1.26) 

x-»® 

lim  x_1R(x)  -  2 ,  (1.27) 

x"»® 

and 

1.485  <c<  1.50  (1.28) 

Since  it  may  be  easily  shown  that  l/  (x)  >  0  from  (1.5),  Abelian 
and  Tauberian  theorems  (Widder  (1941),  pp.  182,  192)  yield  from  (1.24), 
(1.25)  that  (1.26)  holds.  The  equation  (1.27)  is  lnmediate  from  (1.19). 

To  show  that  c  <  1.5  in  (1.28),  taking  the  Laplace  transform  of 
the  second  equation  of  (1.13),  where  I^(x)  =  1.5x  +  1.5  and  repeating 
steps  as  in  (1.21)  -  (1.25)  one  obtains  that 


lim  x-1I^(x)  -  | 
x  *•  • 


(1.29) 


^ jYV^  d“  r  ml  +  +  (e-.« .^1 


t  -2Y  ^  -2J 

4e  '  i  e  v 


f“e^ 
J  .  u 


I  — —  +  r- 

L  v  -  2  3 

2v  s 


V 


By  evaluating  the  contribution  to  the  integral  in  (1.29)  near  0, 
and  comparing  it  to  the  contribution  near  0  to  the  Integral  in  (1.25), 
one  obtains  that  these  are  not  equal. 

To  establish  the  left  inequality  in  (1.28),  a  similar  method  to 
that  used  above  suffices.  This  is  outlined  as  follows.  It  is  easily 
established  that  there  exists  an  a  such  that 


r  1.485x  +  .485  ,  x  >  ■? 

*(x)  =  \  x  1 

satisfies,  for  x  >  j,  ,  x  <  2 


A(x)  ■l  +  x  +  xi  2  ^u>du’ 


(1.30) 


(1.31) 


Taking  Laplace  transforms  of  (1.31)  for  0  <  s  <  1, 


k(s)  =  [  e  8uJt(u)du, 
1 
2 


(1.32) 


one  obtains  an  explicit  solution  for  k(s),  and  then  that  for  0  <  s  «  1, 


k<*>  -§2+§f1 


(1.33) 


where  $  is  expressed  as  a  definite  Integral  which  is  seen  to  be  not 
equal  to  that  in  (1.25).  This  completes  the  proof. 
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II.  Second  Moments;  Covering  Models. 

Let 

M^x)  -  EX2  (x) 
1*2  (x)  =  EY2  (x) 


(2.1) 


Theorem  II. 


lim  x  ^Var  X(x) 
x"*  • 


a  >  0 


where  a  is  obtained  explicitly. 
For  x  >  0, 


Var  Y(x)  -  |(x+l) 


(2.2) 


(2.3) 


Proof .  Squaring  (1.4)  and  taking  expectations,  it  follows  that 

2 

M-(x)  »  2M(x)  -  1  +  -  J  M(t)M(x-l-t)dt  (2.4) 

x  0 

+  I 

2  f  *  ”1 

^(x)  -  2R(x)  -  1  +  J  R(t)R(x-l-t)dt 

•  x 


+  JTT 


with  respective  initial  conditions 
(i)  !^(x)  -  1,  0  <  x  <  i 

(ii)  R^O)  -  1. 


(2.5) 


(2.6) (i) 

(2.6)  (ii) 


Since  R(x)  m  2x  +  1 ,  (2.5)  has  immediate  solution 
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and  hence 


^(x)  -  4x2  +  yx  + 


Var  Y  (x)  «  R^x)  -  R2(x)  -  |(x+l), 


which  is  (2.3). 

To  prove  (2.2),  define 

i(x) 


1, 

cx  +  c 


0  <  x  <  2* 

1,  x  >  i. 


Then,  defining,  for  x  >  1, 


D(x)  =  E (X(x)  -  J5(x)\ 


it  is  obtained  that 


2  f  *  ~2 


D(x)  =  f(x)  +  j  j  (M(t)  -  £(t)) (M(x-l-t)-£(x-l-t))dt 


2  rx“2. 


+  ~j  D(t)dt , 


where 


lim  x”*(f(x)  ■  0, 

x  -*  OS 


and  f(x)  is  expressible  in  terms  of  M(x) . 
Hence  (2.11)  may  be  written 


JL 

D(x)  -  g(x)  +  |  J  2D(t)dt, 
0 


(2.7) 


(2.8) 


(2.9) 


(2.10) 


(2.11) 


(2.12) 


(2.13) 
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where,  by  (1.24),  (2.12) 

lim  x  ^g(x)  ■  b  4  0.  (2.14) 

x  "*® 

Solving  for  the  Laplace  transform  of  D(x)  as  in  section  I,  one  may 
obtain  that  the  asymptotic  behavior  of  D(x)  may  be  established  as 

lim  x  ^D(x)  =  a  >  0.  (2.15) 

x  "*  ® 

Noting  that 

|Var  X(x)  -  D(x)|  -  ^(x)  -  ^(x)  -  (^ (x)-2 2(x)M (x)+*2) 

=  (M(x)-£(x))2,  (2.16) 

it  follows  from  (1.24)  that,  as  x  *•  ®,  (2.17) 

x  * |Var  X(x)-D(x)|  =  o(x),  establishing  (2.2). 

Ill .  Abacus  model. 

Discrete  covering  model  versions  of  those  in  sections  I,  II  are 
given.  A  line  segment  of  length  a  +  1,  with  a  >  0  an  integer,  is  placed 
on  the  x-axis  so  that  the  left  endpoint  is  0.  For  d  >  0  an  integer,  a 
2d-length  "car"  is  set  ("parked")  at  random  upon  the  segment,  such  that 
the  midpoint  of  the  "car"  is  placed  upon  any  of  the  points  (0,1),  (0,2), 
(0,a)  uniformly  with  probability  — .  Independent  of  the  first  "car,"  a 
second  2d- length  car  is  placed  with  midpoint  chosen  uniformly  at  random 
(if  possible)  at  (0,1),  (0,2) , . . . (0 ,a) ,  and  "parked"  (kept  in  that  place) 
if  and  only  if  the  center  of  the  car  is  on  an  uncovered  point  (0,k), 
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1  <  k  <  a  of  the  original  segment,  otherwise  it  is  discarded,  and  a  new 
2d-length  car  is  placed  l.I.D.  as  the  previous  one.  The  process  con¬ 
tinues  until  for  the  first  time  the  points  (0,k)  for  1  <  k  <  a  are  all 
covered.  This  is  a  Renyi  model  of  covering  the  one-dimensional  abacus. 

A  Solomon  model  of  covering  is  similar,  except  that  a  newly  placed 
d-length  car  with  left  endpoint  on  an  integer  point  is  kept  ("parked") 
if  and  only  if  some  portion  of  the  car  covers  a  previously  uncovered 
integer  point  (0,k),  1  <  k  <  a. 

Let,  for  a  >  0  an  integer,  d  >  0  an  integer 
X(a)  ■=  minimum  number  of  2d-length  cars  needed  to  cover  the  [0,a+l] 
interval  in  the  Renyi  model  (3.1) 

Let  Y(a)  »  minimum  number  of  d-length  cars  needed  to  cover 
[0,a+l]  in  the  Solomon  model. 

Denote 

m(a)  =  E(X(a)). 
r(a)  =  E(Y(a) ) . 

<?  (a)  s  E(X(a)-m(a))2 
i*(a)  s  E(Y(a)-r(a))2. 

Theorem  III. 

lim  a  *ra(a)  ■  a.  >  0  (3.4) (i) 

a  1 

lim  a  * r  (a)  ■  a,  >  0 

«  mm  m  “ 


(3.2) 

(3.3)  (i) 

(3.3) (ii) 

(3.3)  (iii) 

(3.3)  (iv) 


(3,4)(ii) 
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lim  a-1a2(a)  -  8  >  0  (3.4)(iii) 

a  -  ®  1 

lim  a-1i?(a)  »  B,  >  0.  (3.4)(iv) 

a-*® 

Proof.  Again  taking  expectations  of  (1.4),  regarding  the  quantities 
there  as  appropriate  for  either  X(a)  or  Y(a),  respectively,  yields  for 
the  Renyi  model, 


,  a-d-1 

m(a)  ■  1  +  —  E  m(A) , 
*  4-1 


(3.5) 


and  m(l)-m(2)-. . .  ,m(d)=m(d+l)«l . 
For  the  Solomon  model. 


a-1 


r(a)  -  1  +  -f-r  E  r(  A) , 


with  r(l)  -  1. 


As  in  section  1,  define,  for  0  <  s  <  1, 


(3.6) 


£ 

K(s)  -  E  e”8  m(  A) . 
A-d+2 


(3.7) 


Since  m(A)  <  A,  K(s)  is  well  defined. 

Clearing  of  fractions  in  (3.5)  and  applying  (3.7)  yields 


*'(•)  + 


i-e'8  \a-e’8)  /  1_e" 


-s(d+2) 


e-»(d+2) 

~  -8.2 

(1-e  ) 


( (d+l)e's-d] 


(3.8) (i) 


,  -s(d+l)  -s(d+2) 

K'(s)  +  il - - -  K(s)  -  - - — =•  [  (d+l)e_8-d]  .  (3.8)  (ii) 

1-e"8  (1-e"8)2 


This  has  solution 


K(s)  ■  exp  2 


■u(d+l) 


1-e 


-u 


dU]| 


exp[-2j 


r»e-u(d+l) 


1-e 


du  j 


Ti 


e-v(^2){(d+i)e-v.dr 

L,,  -v.2  J  dv' 

(1-e  ) 


(3.9) 


For  0  <  s  «  1, 


K(s)  - 


,  2 (d+l)b-l 
s 


(3.10) 


where 


b 


r  r®e'u(d+1)  i  re'v(d+2)  r  -v 

j  expl-2  - — —  du  +  26j  ^ - I\TT  £(**+l)e  <*v, 

"0  lv  l-e‘u  (1-e  ) 


(3.11) 


with 
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,1^  e-u(d+l) 
'  0  ^  (l-e"U) 


r®e-u(d+l) 


—  "y  du  +  I 
u  '  J1  (l-e”u) 


du 


(3.12) 


It  is  noted  that  the  symbols  er^6^y^,  etc.  are  used  to  denote  different 
quantities  in  different  expressions. 

A  similar  result  holds  for  the  Solomon  model,  by  the  same  method. 


outlined  below. 


2 

in  (3.5),  (3.6)  respectively,  and  equating  coefficients  of  a  ,  a,  1, 
respectively,  after  clearing  the  denominators,  one  obtains  for  the 
Renyi  model 

P  a-d-1 

a(a1a+P,)  =  a  +  2  £  (o 

11  ~  j^d+3 

+  2  (d+2+  ^), 


l  +  3 


!>] 


(3.20)  (1 ) 


since 


m(d+2) 


1  + 


d+2 


(3.20) (ii) 


and  one  obtains 


3d2+lld+12 
2  (d+2  )  (2d2+5d+3  ) 


-  0l(2d+l)-l, 


and  for  the  Solomon  model,  the  corresponding  result  is 


(3.21) (i) 


“2 


2 

d+2  ’ 


h 


d 

d+2  ' 


(3.21) (ii) 


The  variances  for  the  two  models  are  obtained  by  squaring  (1.24) 
to  obtain,  denoting 


m^(a)  =  E  (X(a))2  (3.22)(i) 

r2(a)  =  E(Y(a))2,  (3.22)  (ii) 


that,  respectively. 
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0  a-2d-2  a-d-1 

■u(a)  -  2m(a)  -  1+-  E  m(jl)m(a-2d-j&-l)  +  -  E  nu(£)  (3.23) 

*  a  j6-1  a  X-l  * 


with 


^(i)  -  1,  1  <  i  <  d,  ^(d+l)  -  1  + 


0  a-d-2  a-1 

r,(a)  -  2r (a) -1+^r  Z  r( jt)r(a-d- A-1)  +  E  r  (l) 

2  a-H  ^  a+d  2 


(3.24) 


with 


V*>  *  l-  r2(2>  ■  1  +  5« 


Denoting  the  linear  functions,  for  x  >  d 


^(x)  -  bx  +  (2b(d  +  l)-i)  ~  m(x) 
£2(x)  -  kx  +  (k(d+2)-l)  -  r(x) 


(3.25)  (1) 

(3.25)  (11) 


corresponding  to  the  means  for  the  Renyi  and  Solomon  cases,  respectively. 


(1)  s2(a)  -  E(X(a)-jfc1(a))2 

(11)  t2(a)  -  E(Y(a)-42(a))2, 


(3.26)  (1) 

(3 .26)  (11) 


From  (3.23),  (3.24)  respectively,  it  follows  that 


s2  (a) 


,a-d-l 

fx(a)  +  ^  2  (m(v)-jt1(v))(o(a-2d-v-l)-i1(a-2d-v-l)) 

^  v-1 


,a-d-l  , 

+  ~  E  s  £(y) 
v«l 


(3.27) 
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t2  (a) 


2  a-d-2 

f2*a^  +  a+d  E  <r(v)'^(v)Hr(a-d-v-l)-jL(a-d-v-l)) 
v»l  * 


2  a_^  2 
+  r£r  2  t  (v) , 

a+d  v-1 


(3.28) 


where  | f ^ (a) 1  <  M,  i«l,2,  for  some  constant  M  <  ®. 

From  the  Laplace  transform  of  (3.27),  (3.28),  It  may  be  concluded 

that 

a_1s2(a)  -  Ct  <  »  (3.29)(i) 

a-1t2(a)  -  C2  <  ®  (3.29) (ii) 

and  hence  that  (3.4) (ill),  (3.4)(iv)  hold.  This  completes  the  theorem. 


IV.  Random  Car  Size 

One  covering  model  by  randomly  chosen  car  lengths  is  as  follows. 

On  a  curb  [0,a]  a  car  of  length  x  is  placed  with  center  chosen  uniformly 

on  [0,a],  and  x  is  chosen  from  a  probability  distribution  with  density 

f(x).  The  car  is  kept  in  place  (sparked)  if  and  only  if  its  center  is  on  TO, a]. 

A  second  car  is  chosen  I.I.D.  as  the  first,  and  parked  if  and  only  if 

its  center  is  on  an  uncovered  portion  of  [0,a],  The  process  continues 

until  no  further  portion  of  the  curb  is  uncovered. 

Let 

•  x 

F(x)  ■  j  f(y)dy,  with  F(c)"0,  F(c+)  >  0,  for  some  c  >  0  (4.1) 

c 

denote  the  distribution  function  of  one-half  of  the  car  length,  and 

M(a)  ■  mean  total  number  of  cars  required  to  first  cover  [0,a].  (4.2) 
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Then  by  considerations  as  in  (1.4)  for  a  >  c. 


2 ' a-u 

M(a)  «  1  +  -J  f (u)du  M(y)dy  (4.3) 

a  0  0 

and  M(x)  *  1,  0  <  x  <  c. 

A  formal  solution  for  the  Laplace  transform  of  may  now  be  obtained. 
It  is  assumed  that  M(a)  is  finite.  This  can  be  guaranteed,  e.g., 
by  setting  F(c)  «  0  for  some  c  >  0. 

Assume  that 


2 

y  f (y)dy  <  ».  (4.4) 

0 


Denote 


*  00 
'  0 


$(s)  •*  I  e"8yf(y)dy. 


(4.5) 


L(s)  =  e  SyM(y)dy. 


(4.6) 


Then  multiplying  (4.3)  by  e"sa  and  integrating  from  0  to  ®  yields 


+  m*i 


L'(s)  +  L(s) 

s 


e~8Cq+sc)  +  28(s)(l-e~SC) 

2 

s 


(4.7) 


with  solution 


w 

!•(.)-. xpfj  '“iisl  du  .  Lj\2  du+26  f  e-vctl-Wc)4-2»(v)(l-e-vc)  - 

-  _  u  J  i  r  „  u  _  2  J 

8  V  tr 


dv 


(4.8) 
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with 


6 


,1 

,l 

0 


i  frO-i 


da  + 


du. 


For  0  <  s  «  1, 

«•)=  -T  +  Ui^i 

8  8 


(4.9) 


(4.10) 


where  here 


V 


J  exp[-2j  £^du+26_|j 

0  v 


e-vc(14vC)+2_|.(v)^e:VC)1  dv 

v 


(4.11) 


and 

CO 

n  -  '  yf(y)dy.  (4.12) 

vo 

Theorem  IV.  Under  the  hypotheses  of  section  IV, 

lim  a_1M(a)  -  y.  (4.13) 

a  -*ob 

Proof.  This  follows  from  (4. 10) and  a  Tauberlan  theorem  (Widder 
1941),  pp.  182,  192),  since  M(a)  is  of  bounded  variation. 

Remark  I.  A  simple  approximation  to  y  may  be  obtained  by  finding 

the  approximate  linear  solution  to  (4.3).  Assuming  a  solution  to  (4.3) 

2 

of  form  aa  +  8,  one  equates  coefficients  of  a  ,a,l  for  a  large,  in 

a  r--a-u  i 

a(aa48)  -  a  +  2j  f(u)du)  j  (ay+0)dy  +  c  j  ,  (4.14) 

0  c 
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that  is,  solve  for  or, 9  in 

cra^  +  9a  *’  a  +a[a^  -  2ap  +  -  c^]  +  29[a-u-c]  +  2c 


(4.15) 


where 


to  obtain 


T  2 

v>2  *  J  y  f(y>dy. 


2u  +  4c 


2  2 
c  +  4(1  +  4 pc  -  n 


,  9  ■  2cfli  -  1 


(4.16) 


(4.17) 


The  properties  of  this  approximation  remain  to  be  investigated.  It  is 

conjectured  to  be  a  good  upper  bound.  When  c  *  |i  *  j,  “  £  * 

3 

corresponding  to  unit  length  cars,  a  m  which  is  a  good  upper  bound 
in  the  Renyi  covering  model  of  section  I. 


Remark  II.  Let 

W(a)  ■  total  number  of  random  size  cars  required  to  cover  [0,a],  (4.18) 
Let 

M^a)  -  EWZ(a)  (4.19) 


By  the  usual  argument  as  in  (1.4), 

„  ?a  ,-a-u 


KL(a)  -  2M(a)  -  l+|j  f(u)dUj  M(y)M(a-u-y)dy 
i  0  0 


2 ?a  c a-u 

+  -\  f  (u)du  J  M,(y)dy 
a‘  0  0 


(4.20) 


and  Mj(y)  -  1,  0  <  y  <  c. 


Since  this  expression  is  similar  in  form  to  the  second  moment 
expressions  for  the  models  of  previous  sections,  it  is  plausible  that 


VarfW(a)]  -  Mj (a)-(M(a))2  -  E  (W(a)-A(a))2  (4.21) 

where 

1(a)  -  ya  +  2nY-  1, 

and  hence  it  is  conjectured  that,  under  mild  moment  conditions, 

lim  a-1Var[W(a)]  -  T\>  0, 

a  -*  ® 

where  T|  is  a  finite  constant. 

V.  Alternating  car  size 

On  a  curb  [0,x] ,  a  car  of  length  1  is  parked  in  accord  with  a 

Renyi  model  (i.e.  the  center  must  be  on  [0,x])  or  a  Solomon  model 

(i.e.  some  portion  must  be  on  [0,x]).  The  next  car  candidate  in  the  same 
model  is  of  length  a  >  1,  and  in  the  Renyi  model  is  parked  if  and  only  if, 
the  center  falls  on  an  uncovered  portion  of  [0,x] ,  and,  for  a  >  2,  the 
car  does  not  completely  contain  the  car  of  length  1.  In  the  Solomon 
model,  some  portion  of  the  a-length  car  must  fall  on  an  uncovered  portion 
of  [0,x],  and  also  must  not  completely  contain  a  car  of  length  1, 
otherwise  it  is  discarded.  If  an  a-length  car  is  discarded  in  either 

model,  another  a-length  car  is  selected  I.I.D.  as  the  prior  one, 

until  one  is  parked,  or  the  Interval  [0,x]  is  completely  covered,  at 
which  time  the  process  stops.  If  the  a-length  car  is  parked,  and  some 


(4.22) 


(4.23) 


i 

i 
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portion  o£  [0,x]  is  uncovered  in  either  model,  the  next  car  to  be 
parked  is  one  of  length  1,  followed  by  one  of  length  a,  and  so  on, 
alternately,  until  [0,x]  is  completely  covered  in  either  model,  at 
which  time  the  process  stops. 

Let 

X^(x)  -  minimum  total  number  of  alternate  size  cars  required  to 
cover  [0,xj  in  the  Renyi  model  starting  with  a  car  of  length  1,  and 
X^(x)  =*>  same  quantity  starting  with  a  car  of  length  a. 

X ' " ^ (u)  -  minimum  total  number  of  alternate  size  cars  required  to 
cover  [0,u]  in  the  Renyi  model  when  the  first  car  to  be  parked  was  an 

a-length  car,  parked  at  [u,u+a]  on  a  total  interval  TO,x],  for  x  >  u  +  a. 

(2) 

X  '(u)  =  minimum  total  number  of  alternate  size  cars  required  to 
cover  [0,u]  in  the  Renyi  model  when  the  first  car  to  be  parked  was  a 
1-length  car,  parked  at  [u,u  +  l]  on  a  total  Interval  [0,x] ,  for 
x  >  u  + 1. 

Let  (5.2) 

Y^(x),  Y^x),  Y^\u),  Y^(u)  be  the  corresponding  quantities  for 
the  Solomon  model. 

Let  (5.3) 

r^x)  -  E(Yt(x))  i-1,2 
r(i)(x)  -  E(Y(i)(x))  i-1,2. 

By  considerations  as  in  the  previous  sections,  by  expressions 


similar  to  (1.4),  it  follows,  for  the  Solomon  model,  for  example. 
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2  r° 

r2(x)  "  7+a  J  1+r1(x-“-a)du 


+  J  du[l  +  (^2)(rl(u)+r<2)(x-u-a))].  (5.4) 


rl<x)  »  du(l+r2(x+l-2a-u) 


J+T  £“[*  +(M^f!)Ct2("+2-2*)  +r(I>(«-»-l)). 


1^(0)  »  r2(0)  ■  1. 


Also,  an  approximation  which  may  be  good  for  large  x  is 


(5.4) (11) 


n* ,««.  j.”  («»).,« 


*rioaa->««. 


(5.4) (ill) 


Similar  expressions  hold  for  the  means  in  the  Renyl  models,  with 

slightly  different  initial  conditions,  and  the  cases  1  <  a  <  2  and  a  >  2 

are  treated  separately.  It  may  be  possible  to  obtain  linear  solutions  to 

(5.4)  for  x  large,  but  the  computation  is  tedious,  and  of  unknown  merit. 

The  exact  solution  of  (5.4)  by  Laplace  transforms  appears  possible  by 

power  series  or  a  Volterra-like  expansion,  and  an  asymptotic  solution 
-2 

with  leading  term  s  is  sought,  corresponding  to  asymptotic  linearity, 
since  it  is  conjectured  that  r^x),  i-1,2,  grow,  for  x  —  • 


r£(x)  ~  OjX  +  $1,  i-1,2. 


(5.5) 


This  is  because  r  (x)  <  r  (x)  <  r(x),  where  r(x)  is  the  mean  number  of 
«  I 

unit  length  cars  required  to  cover  [0,x],  and  rfi(x)  is  the  mean  number 
of  a- length  cars  required  to  cover  [0,x]  in  the  Solomon  model. 


VI .  Abacus  Packing. 


A  one -dimensional  abacus  packing  (s  parking)  model  on  a  discrete 
grid  is  considered.  A  car  of  length  d,  where  d  is  an  odd,  positive 
integer,  is  parked  horizontally  with  center  chosen  uniformly  from 
YpSO )J  f^,0),..V-  where  a  is  a  positive  integer,  a  >  d. 

A  second  car  of  length  d  is  chosen  I.I.D.  as  the  first,  and  parked  if 
and  only  if  it  does  not  intersect  the  first  car,  and  its  end  closest  to 
the  first  car  is  at  least  1  unit  away,  otherwise  it  is  discarded,  and 
another  car  chosen  I.I.D.  as  the  first,  and  so  on,  until  it  is  parked 
or  no  further  cars  may  fit  in  the  prescribed  manner. 

Let  (6.1) 

X(a)  ■  total  number  of  d-length  cars  (d  an  odd,  positive  integer) 
which  may  be  parked  on  {(£,0)  1  <  t  <  a},  with  center  at  (k,0),  with  k 
chosen  uniformly  from:  ^ -  <  k  <  a  -  (^jp). 

Then,  conditional  on  the  center  k  chosen  as  above. 


x(a)  -  { 


1+  X’(k-(~))  +  X"(a-k-(~))  if  <k  <a 


(±tL) 

*  2  ' 


1  +  X’(a-d-l)  if  k  -  y  or  k 


a-(^) 


(6.2) 


where  X’(r),  X"(r)  are  I.I.D.  as  X(r). 
Let,  for  r  >  0  an  integer. 


m(r)  s  EX(r) 


(6.3) 
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Then,  from  (6.2),  on  taking  expectations,  one  may  obtain 


2 

m(a+d-l)  -  1  +  -  E  m(i)  (6.4) 

a  i*  1 

with  m(f)  =  1 ,  d  <  1  <  2d,  m(£)  “0,  £  <  d. 

To  obtain  the  asymptotic  mean  number  of  d-length  cars  which  may  be 
parked,  one  proceeds  as  before  using  the  Laplace  transform. 

Define 

“  i 

K(s)  =  Em(£)e"S  (6.5) 

fed 

Then  (6.4)  yields  that,  summing  from  J&-1  to  °°, 


,  -s(d+l)  N  -e'sd 

K'(s)+'d-l+— - - -  J  K(s)  »  - — r  (6.6) 

l-e~S  ^  (l-e-S; 


with  solution 


r  '’a,e"uCd+l)  1  f  r  f  e-u(d+l)  ip  e“dv  1 

K(s)*exp^-(d-l)s+2  -  .7.^  du  !  expj_(d-l)v-2j  - — —  du j  j  - ^TTJ  dv’ 

s  (1-e  )  '  J  v  1-e  *  (!  ) 


(6.7) 


so  that  for  0  <  s  «  1, 


K(s) 


*  c  ,  c(d+3)-l 


2  s 


(6.8) 


where 


P®  r  r®  e-u(d+l)  ip  e-dv 

c  -  I  exp|_(d-l)v-2J  —  du  +  2flj|_  _v^  J 


v  (1-e  “) 


/i  -v»2 
(1-e  ) 


dv 


(6.9) 
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with 


r»e-u(d+l)  plr  e~u(d+l)  1 

I  -  du  +  I  — - —  | 

Jl  (l-e"U)  0  L  (i-e“U)  u  J 


du. 


(6.10) 


A  Tauber ian  argument  (Widder  (1941),  pp.  182,  192),  yields, 
for  j l  large 


m(f)  ~  cf  +  c(d +  3)  -  1. 


(6.11) 


Remark  1.  An  approximation  to  m(f)  in  (6.11)  may  be  obtained  by 
obtaining  the  linear  solution  to  (6.4),  as  follows. 

Writing  (6.4),  for  a  >  2d  +  3  as 

r  a-2  i 

a  m(a+d  -  1)  -  a  +  2  I  E  m(f)|  +  2(d+l)  (6.12) 

u  i-2d+l  J 

and  trying  a  solution  of  the  form 

r(f)  -  al  +  3  (6.13) 

2 

and  equating  coefficients  of  a  ,  a,  1  respectively  yields  the  solution 

,  31  +  3  -  d 

“  4d  +  3  '  (6.14) 

Remark  2.  The  properties  of  a  linear  solution  depend  on  d. 

The  case  d  »  1  is  given,  as  in  Weiner  (1980). 

The  equation  (6.4)  becomes 

2  a-^ 

m(a)  «  1  +  -  E  m (1)  (6.15) 

*  1.1 


^ - - .  . 
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with  m(l)  -  m(2)  -  1. 


Lemma  2 .  If  for  i  >  r  +  2 ,  for  some  r  >  1 , 


A(i)  -  \  E  A(k) 
k-r 


(6.16) 


and  A( f.Q+1 )  >  0,  A(fQ+2)  >  0,  £Q+1  >  r  +  2,  then  A(£)  >0,  l  >  +  1. 


Proof. 


A(  Xg+2 )  >  0  implies  Z  A(k)  >  0, 

k-r 


(6.17) 


so  that 


A<V3)  “  T+T  L  ^  &<k>  +  A(Jtn  +  l)]  >  0. 

0  k-r 


(6.18) 


Now  A(£q  +  2)  >0,  A(jJq  +  3)  >  0,  and  a  repeat  of  the  argument 
yields  a  proof  by  induction. 


Lemma  3.  For  l  >  5, 


3£  +  2  .  Ul+3 

-  <  m(£)  <  - 

7  “9 


(6.19) 


Proof.  Eq.  (6.15)  may  be  written 


2  2  W 

m(£)  -  1  +  7  +  7  E  m(k), 
1  1  k-2 


(6.20) 


For  A  >  4,  it  is  imoediately  verified  that  the  function 


...  _  4A  +  3 
r(f)  =  — g — 


(6.21) 


satisfies 


2  2 

1  +  7  +  7  E  r  (k) . 
£  X  k-2 


(6.22) 


Denoting 


Aj(£)  -  r ( A)  -  m( A) , 


(6.23) 


Akf2  Vk> 


(6.24) (i) 


and  a  computation  yields  that 


^(5)  >  0,  ^(6)  >  0. 


(6.24) (ii) 


Hence  lemma  1  applied  to  (6.22)-(6.24)  yields  the  right  side 
inequality  of  (6.19). 

For  A  >  5,  (6.15)  may  be  written 


4  2  4 

m( A)  ■1+7  +  7  E  m(k) 
1  1  k-3 


(6.25) 


and  for  A  >  5,  the  function 


V(A)  s  H_+l 


(6.26) 


satisfies 


4  2 

V(A)"  1  +  7  +  7  2  V(k). 
A  l  k-3 


(6.27) 
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Denoting 

Z^(i)  -  m(i)  -  V(je),  (6.28) 

then 

2  t-2 

^2<je)  "  ^  ^2  (k>  (6.29) 

and  a  computation  yields  that 

^(5)  >  0,  ^(6)  >  0,  (6.30) 

so  that  lemma  1  applied  to  (6.28)-(6.30)  yields  the  le£t  inequality  of 
(6.19). 

By  arguments  for  the  second  moments  as  in  the  previous  sections, 
it  may  be  seen  that  the 

lim  x-1Var  X(x)  *»  y  >  0  (6.31) 

x-® 

VII .  Remarks . 

Remark  1.  Using  Weiner  (1979),  it  is  easily  seen  from  higher 
asymptotic  moment  computations,  obtained  by  taking  expectations  of 
higher  powers  of  (1.24)  that  a  central  limit  theorem  holds  for  each 
model  in  sections  I,  II.  One  example  is  the  following. 

Let  (7.1) 

X(x)  ■  total  number  of  cars  required  to  cover  [0,x]  in  either  a 
Renyi  or  Solomon  model  of  sections  I,  II. 


Denote 


lim  x-1E(X(x))  M>0 
x  "• 00 


<7.2)(i) 


(7.2) (ii) 


lim  x“1Var(X(x))-  o'  > 

x-*«9 

Theorem  V.  For  a  <  b, 

lim  l[.  <  g  <  bl  .  -L 

x“*®  aVx’  V2tt 


0. 


(7.3) 


Remark  2. 

The  same  considerations  are  conjectured  to  suffice  to  establish  a 
central  limit  theorem  for  models  in  III,  IV,  VI,  and  a  bivariate  central 
limit  theorem  for  the  alternate  car  size  model  of  section  V. 

Remark  3.  It  is  conjectured  that  a  two  dimensional  "Palasti" 
result  holds  for  the  covering  models  of  sections  I,  II.  That  is,  for 
unit  square  cars  dropped  uniformly  on  a  rectangular  parking  area  of  size 
x  X  y,  let  (7.4) 

M(x,y)  *  mean  total  number  of  unit  square  cars  required  to  cover 
x  X  y  under  a  Renyi  or  Solomon  model. 

The  Palasti  conjecture  for  this  model  is  then 

lim  (xy)_1M(x,y)  *  .  (7.5) 

x,y  "•  08 

See  Weiner  (1978)  for  a  heuristic  discussion  of  packing  models  in  the 
plane. 
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